AC susceptibility investigation of vortex dynamics 
in nearly-optimally doped REFeAsOi^F,,, superconductors (RE = La, Ce, Sm) 
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Ac susceptibility and static magnetization measurements were performed in the nearly-optimally 
doped LaFeAsOo.gFo.i and CeFeAsO0.92F0.08 superconductors, complementing earlier results on 
SmFeAsOo.8F .2 [Phys. Rev. B 83, 174514 (2011)]. The magnetic field - temperature phase 
diagram of the mixed superconducting state is drawn for the three materials, displaying a sizeable 
reduction of the liquid phase upon increasing T c in the range of applied fields (H < 5 T). This result 
indicates that SmFeAsOo.sFo.2 is the most interesting compound among the investigated ones in 
view of possible applications. The field-dependence of the intra-grain depinning energy Uo exhibits 
a common trend for all the samples with a typical crossover field value (2500 Oe < H cr < 5000 
Oe) separating regions where single and collective depinning processes are at work. Analysis of the 
data in terms of a simple two-fluid picture for slightly anisotropic materials allows to estimate the 
zero-temperature penetration depth A a f>(0) and the anisotropy parameter 7 for the three materials. 
Finally, a sizeable suppression of the superfluid density is deduced in a s two-gap scenario. 

PACS numbers: 74.25.Uv, 74.25.Wx, 74.70.Xa 



I. INTRODUCTION 

Almost four years after the discovery of high- 
temperature superconductivity in Fe-based pnictidcsj 1 " 
several questions are still open both on fundamental as- 
pects and on possible technological applications. No clear 
and exhaustive explanations for the precise pairing mech- 
anism have been given yet. Incontrovertible features of 
isotopic effects cannot entirely rule out a partial role of 
the lattice in the coupling process^ but it is essential 
to stress that T c 's values are indeed too high to allow 
only a conventional pairing to be at worki^ Strong evi- 
dences for a multi-band s scenario similar to that pro- 
posed for magnesium di-boride^ have been reported, to- 
gether with theoretical support from a number of tech- 
niques in materials belonging to the 1111 family^— to 
122 compounds^ and to 11 chalcogenidesJ^ 

At the same time, it is still not clear whether these 
novel compounds will be helpful as valid technological 
tools. The answer to these problems will essentially come 
from specific measurements of critical current densities in 
relation to the features of the grain boundaries. Detailed 
studies of the practical applicability of granular supercon- 
ductors, anyway, cannot leave aside precise investigations 
of the fundamental intrinsic properties of the materials. 
In this respect, the determination of the so-called irre- 
versibility line is of the utmost importance. Such line, in 
fact, delimits the region in the magnetic field - tempera- 
ture phase diagram where the dissipationlcss property of 
the thermodynamical superconducting phase is preserved 
even in the presence of a partial penetration of magnetic 
field inside the material. In a previous work on a powder 
sample of optimally-doped SmFeAsOo.sFo.2 (see Ref. [l3l ) 
the determination of the intrinsic irreversibility line was 
performed by means of the ac susceptibility technique. 



Parameters like the superconducting critical temperature 
T c are anyway well known to be strongly dependent on 
the considered RE ionJ^ The question of whether (and 
how) the properties of the phase diagram of the flux lines 
and of the relative irreversibility lines too are indeed de- 
pendent on the RE ion is then of extreme importance. 



In this paper we report on the phase diagram of the 
flux lines in three powder samples of REFcAsOi-^Fa; su- 
perconductors (RE = La, Ce, Sm) under conditions of 
nearly-optimal doping. In particular, much attention is 
devoted to the determination of the irreversibility line 
and its dependence on the RE ion. The features of pin- 
ning mechanisms and, in particular, of the characteris- 
tic depinning energy barriers are investigated in detail 
within a thermally-activated framework. Results are in- 
terpreted and analyzed by distinguishing two magnetic 
field regimes characterized by single and collective depin- 
ning processes. Together with a simple two-fluid model 
and a proper normalization, this allows us to make ex- 
perimental data collapse on the same curve, a feature 
indicative of a common underlying mechanism indepen- 
dent on the precise material. Reliable estimates of the 
zero-temperature penetration depth A a {,(0) and of the de- 
pendence of the anisotropy parameter 7 on the RE ion 
are finally given, together with some interpretations of 
the observed behaviour in terms of a two-band s model. 
Data relative to the Sm-based sample have already been 
presented in Ref. [HI and will be reported also here for 
the sake of clarity and completeness. 
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II. EXPERIMENTALS. ASPECTS OF DC 
MAGNETIZATION AND AC SUSCEPTIBILITY 

Powder samples of REFeAsOi-^F^. (with RE = La, 
Ce, Sm and nominal F~ contents x = 0.1, 0.08 and 0.2, 
respectively) were synthesized as described in previous 
works . 13 i 15 ~ — Static magnetization Md c and ac suscepti- 
bility Xac measurements were performed by means of a 
Quantum Design MPMS-XL7 SQUID magnetometer and 
of a MPMS-XL5 SQUID susceptometer, respectively. In 
the latter case a small alternating magnetic field H ac with 
frequency u m is superimposed to a much higher static 
magnetic field H. Measurements were always performed 
in field-cooled (FC) conditions with H ac = 0.0675 — 1.5 
Oe parallel to /io-ff, which varied up to 5 T, while v m 
ranged from 37 to 1488 Hz. 

Mdc/ H vs. temperature (T) raw data obtained un- 
der FC conditions at H = 5 Oe are shown in Fig. m 
For the sake of clarity the curves have been reported af- 
ter subtracting slight spurious contributions, leaving data 
above the superconducting onset at a constant zero offset 
value. T c (0) is defined as the critical temperature T c for 
H — > Oe. Such values are obtained as the zero intercept 
of a linear extrapolation of data below the diamagnetic 
onseti 3 - and reported in Fig. [1] in correspondence of the 
diamagnetic onsets of the three samples. The saturation 
absolute value of the diamagnetic signals can be evalu- 
ated around (0.27 ± 0.05) in 1/Att units. These values are 
typical for fully-superconducting powder samples in mea- 
surements under FC conditions, the observed reduction 
originating from the sample's geometrical properties and 
morphology. In the case of the CeFeAsO0.92F0.08 sam- 
ple, a sizeable paramagnetic contribution from the Ce 
sublattice can be clearly discerned already at such low 
value of magnetic field, mainly due to the high value of 
the Ce 3+ magnetic moment. In fact, the fitting procedure 
to the experimental data at several values of H already 
described in a previous worki^ yields to fj-ce — 2.1 [Ib 
and nsm — 0.3 \xb (raw data not shown). 

Curves reported in Fig. [T]show quite sharp supercon- 
ducting transitions. The slight roundness of the onset 
may be due to several reasons, both intrinsic (for instance 
the influence of superconducting fluctuations near T c , as 
already discussed in Ref. fl5T ) and extrinsic (slight chem- 
ical inhomogeneity of the F~ doping ions, distribution 
of the geometrical size of grains). The signal, moreover, 
comes from differently oriented grains in the powder sam- 
ple and also this fact may contribute to a broadening of 
the transition. The corresponding powder-averaged up- 
per critical field 

(H c2 (T)) pwd ~ ^H c2Mllab (T) + ^H c2 . mc (T) (1) 

was deduced for the three samples by examining the field 
dependence of T c , as reported in Ref. [l]| 

In Fig. [2] some typical Xac vs. T curves 
for CcFcAsO0.g2F0.08 arc displayed (raw data of 
SmFeAsOo.8Fo.2 have been reported previouslyi 3 .) . Raw 
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FIG. 1: (Color online) Mdc/ H vs. T curves (volume units) in 
FC conditions at H = 5 Oe. The estimated T c (0) values are 
indicated in correspondence of the diamagnetic onsets. 



data for LaFeAsOo.gFo.i are qualitatively similar to those 
in Fig. [2] and are not presented. The x'ac vs - T curves 
can be described in terms of a mixed-state shielding re- 
sponse with some degree of distortion occurring in those 
T regions where sizeable contributions to the imaginary 
component appear. x'ac vs. T curves are composed of two 
main peaked contributions: a narrow peak appears just 
below the diamagnetic onset while a much broader one is 
present at low T . This is quite a common phenomenol- 
ogy in superconducting powder samples . 18 > 19 The high-T 
peak is generally associated with the power absorption 
due to losses inside the single grains, while the broad 
low-T peak can be associated with the generation of weak 
Joscphson-likc links among the different grains. 

The three main features shown in Fig. [2j namely the 
diamagnetic onset in x'ac an( i the two peaks in x'aci are 
strongly shifted to lower T on increasing H for all the 
samples. The -ff-dependence of the diamagnetic onset in 
Xac- m particular, is much more marked than in the case 
of the diamagnetic onset in Md c /H vs. T curves. The 
described behaviour can be directly associated to vortex 
dynamics and to the precise features of the irreversibility 
line i 13 ' 20 " — When dealing with an electromagnetic wave 
impinging on a type II superconducting material, in fact, 
one has to carefully take into consideration typical spa- 
tial penetration lengthscales. The general treatment of 
the problem has been considered in several papers^ - — 
In the presence of vortices, in particular, the overall pen- 
etration depth A ac for the radiation can be generally ex- 
pressed as the quadrature sum of Al (representing the 
London penetration depth of the superconductor) and 
Ac (the so-called Campbell penetration depth), namely 



A ac — \/A L + A c . 



(2) 
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FIG. 2: (Color online) Left: Xac vs. T raw data for CeFeAsO0.92F0.08 (H ac = 1.5 Oe and v m = 478 Hz) at different H values. 
The imaginary and the real components are shown in the upper and lower panels, respectively (volume units). Right: zoom of 
the onset region of the curves already reported in the left panels relative to x'ac an d to the first derivative of x'ac with respect 
to T (upper and lower panels, respectively. See text for further details). 



Ac can be generally expressed ai 



,2* 



Ac = 



47ra 



(3) 



where the Labusch parameter a mimics the curvature of 
the potential well associated with the pinning centers in 
a harmonic approximation. In other terms, a quantifies 
the average elastic restoring force density of the pinning 
centers acting on the flux lines (FLs). 

The case of high effectiveness of the pinning mecha- 
nisms (a — > 00) ideally corresponds to a condition where 
FLs arc completely fixed. This condition is, by definition, 
the so-called glassy phase of FLs where typically irre- 
versible processes develop^ By qualitatively considering 
Eqs. ^ and © under these circumstancies, one notices 
that the penetration of the electromagnetic wave is only 
governed by the London penetration depth (A ac — > A^). 
As a result, the electromagnetic wave is shielded by the 
superconductor leading to a diamagnetic response in x'ac- 
On the other hand, in the opposite case of completely in- 
effective pinning the Labusch constant can be considered 
as a vanishing quantity. In this case, the FLs are in the 
so-called liquid state and are substantially free to move 
generating dissipation. The condition a — > yields to 
Ac — > 00 or, at least, to Ac 3> R in the case of a pow- 
der sample (i? is the typical grain size) and no shielding 
can be detected even in the presence of a robust thermo- 
dynamical superconducting phase. As a result, x'ac = 
even if Md c /H already takes negative diamagnetic val- 
ues. The onset of the diamagnetic response of the ma- 
terial in x'ac vs - T curves can then be interpreted as a 
crossover between the two described phases of the FLs 
and its .ff- dependence is a good choice in order to define 
the irreversibility line. 

A slightly different criterion for determining the irre- 



versibility line can be formulated by considering the in- 
trinsic dissipative response inside the grains . 13 ' 22 ! 30 " — In 
particular, by now considering the x'ac vs - ^ curves, one 
can denote by T p the position of the intra-grain maxi- 
mum which is typically found just few-K below the dia- 
magnetic onset in x'ac- This peak can be interpreted as 
arising from a resonating absorption of energy when the 
frequency of the radiation matches the inverse character- 
istic relaxation time l/r c of the vortices in the pinning 
potential dip, namely 



1. 



(4) 



The .ff-dependence of T p can then be chosen in order 
to describe the irreversibility line and in the following 
we will be referring to this as the x'ac~ CT ^ erl0n - Within 
a Dcbye-like relaxation framework, T p almost coincides 
with the characteristic temperature of the correspond- 
ing peak in the derivative of x'ac with respect to T. 
This was experimentally verified in all the three inves- 
tigated samples^ 3 - as explicitly shown only in the case of 
CeFeAsO0.92F0.08 in the right panel of Fig. [21 



III. PHASE DIAGRAMS AND DEPINNING 
ENERGY BARRIERS: MAIN RESULTS 

The FLs phase diagrams for the three examined sam- 
ples are shown in Fig. [3] where both {H C 2(T)) pw d (full 
symbols) and the irreversibility lines (opened symbols) 
are plotted as a function of the reduced temperature 
t = T/T c (0). 



(Hc2{T)) pw d vs. t curves as determined from Md c 
data clearly display linear trends as a function of t. At 
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FIG. 3: (Color online) Phase diagrams for the FLs in the in- 
vestigated samples. Data for SmFeAsOo.sFo.2 have already 
been reported in Ref. [l3|. Open symbols track the irre- 
versibility line as deduced from data taken at u m = 37 Hz 
(dashed-dotted lines are guides for the eye). Full symbols 
track (H c u(T)} pw d (continuous lines are linear fits). 



from the relation 3 ^ (see Tab. HJ) 

(H c 2(0))p W d = 



(6) 



Results from Xac data will now be considered. The 
x" c -criterion was chosen for the determination of the ir- 
reversibility line. Data relative to the derivative of x'ac 
were always analyzed due to their more favourable signal- 
to-noise ratio. A dependence of the T p value on v m was 
detected at all the values of the applied H . For this rea- 
son, data for the lowest accessible value v m = 37 Hz were 
chosen in order to draw the irreversibility lines in Fig. |3] 
It is possible to observe that the extension of the liquid 
phase of the flux lines is progressively reduced by the in- 
crease of T c in the explored H range. This observation 
makes Sm-bascd materials rather interesting in view of 
possible technological applications. 

In order to better compare the behaviour of the three 
samples from a fundamental point of view, the irre- 
versibility lines reported in Fig. [3] are presented again in 
Fig. [4] after normalizing the field values by (H C 2(0)} pw d- 
In all the three samples the irreversibility line can be 
described by means of a power-law function 



low-f values, anyway, a slight upper curvature can be 
detected for all the samples, a feature possibly associ- 
ated with two-band superconductivity^ Data relative 
to CeFeAsO0.92F0.08 are limited to the low- J? values due 
to the dominant paramagnetic contribution arising from 
the Ce 3+ sublattice which fully covers the superconduct- 
ing response for noH > 2 T. The slope of a linear fit 
to the data allows to estimate the upper critical field 
(H C 2(0)) P wd extrapolated at T = K under the simplified 
assumption of single-band superconductivity through the 
Wcrthamer, Hclfand and Hohcnbcrg (WHH) relation^ 



H 



(H c 2(0)) p wd 



0.95- 



T 



-1 3/2 



Tc(0) 



(7) 



characterized by the exponent j3 = 3/2 (see the con- 
tinuous lines in Fig. 0|. This is a typical result in 
high-T c superconductors ) 13 ' 20 ' 37 even if slightly differ- 
ent functional form have been reported, for instance, 
in the case of Ba(Fei_ K Co 2: )2As2 single-crystals.— Here 
the value 0.95 (common for all the samples) phenomeno- 
logically accounts for the discrepancies at low magnetic 



(H C 2(0)) pwd * 0.7 x T c (0) 



d{H c 2)pwd 



dT 



(5) 



T~T c (0) 



An overall correlation is observed between the gradual in- 
crease of T c (0) and a corresponding increase in the slope 
value and, accordingly, in the extrapolated (H C 2(0)) pw d- 
This, in turn, leads to a steady decrease of the extrap- 
olated value of the powder-averaged Ginzburg-Landau 
(GL) coherence length (£{0)) pwd at T — K, calculated 



TABLE I: Critical temperature T c (0) for H -> Oe and cor- 
relation length (£(0)} pw d for T — > K as deduced from Eqs. 
© and @. 



Sample 


T c (0) (K) 


(f(0)) pw 


i (A) 


SmFeAsOo.gFo.2 


52.3 ± 0.25 


11.0 


± 


0.1 


CeFeAsO0.92F0.08 


34.25 ± 0.25 


15.7 


± 


0.3 


LaFeAsOo.gFo.i 


20.2 ± 0.25 


26.3 


± 


0.5 
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FIG. 4: (Color online) Representation of the irreversibility 
lines relative to the investigated samples (already reported 
in the phase diagrams, see Fig. [3} where Hi rr (T) has been 
normalized with respect to the relative value of (H C 2(0)) pw d- 
Continuous lines are best-fits according to Eq. 0. 
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fields when defining the irreversibility line from the x'a C - 
criterion. It should be noticed that, when plotted on this 
different scale for the different samples, the extension of 
the liquid phase as a function of the RE ion shows a 
trend opposite to what displayed in Fig. [3] This inter- 
esting feature will be recalled and discussed later in Sect. 

El 

By now focussing on the ^-dependence of T p 
one can notice that, similarly to what observed in 
SmFeAsOo.8Fo.2r^ the quantity 1/T p displays a logarith- 
mic dependence on v m . This is clearly shown in the inset 
of Fig. EJfor CeFcAsO0.92F0.08 (H ac = 1.5 Oc, = 1.5 
T) even if the phenomenology is well verified for all the 
samples at all the H values. In particular, data can al- 
ways be fitted within a thermally-activated framework 
by the expression 



T p {v m ) (U {H)) pwd \ u J 

(see the fitting function in the inset of Fig. HJ. One 
can recognize that the logarithmic behaviour of 1/T p is 
mainly controlled by the powder-averaged fitting parame- 
ter {Uo(H)} pw d 7 playing the role of an effective depinning 
energy barrier in a thermally-activated flux creep model. 
The parameter i/q in Eq. (JSJ represents an intra-valley 
characteristic frequency associated with the motion of the 
vortices around their equilibrium position in the pinning 
centers. 




0.01 0.1 1 10 



H„H (T) 

FIG. 5: (Color online) Main panel: H-dependence of the de- 
pinning energy barriers (Uo(H)) pw d in the three investigated 
samples. Data relative to SmFeAsOo.8Fo.2 have already been 
reported in a previous work.— Continuous lines are best-fits to 
data according to a simple 1/H dependence. Inset: logarith- 
mic ^-dependence of 1/T P in CeFeAsO0.92F0.08 (H ac = 1.5 
Oe while ^o^ = 1.5 T). The continuous line is a best-fit func- 
tion according to Eq. (JSJ. 

The results of the fitting procedure to T p (y m ) data 
according to Eq. (JSJ have been reported in the main 
panel of Fig. [5] for all the samples. A strongly-marked 



if-dependence of (Uo(H)) pw d can clearly be discerned. 
Remarkably, a common trend is displayed for all the 
samples, making it possible to guess a common under- 
lying mechanism. Beyond an overall sizeable reduc- 
tion of (Ua(H)) pw d with increasing H, a sharp crossover 
(at field values 2500 Oe < H cr < 5000 Oe common 
to all the samples) between two qualitatively different 
(Uo(H)) pw d vs. H regimes is observed for each sample. 
At low fields H < H cr the depinning energy (Uo(H)) pw d 
is found to be only slightly dependent on H, while 
for values H > H cr a trend (Uo(H)) pw d oc 1/H can 
be discerned. The 1/H regime for depinning energies 
(Uo(H)) pw d was indeed observed in different supercon- 
ducting materials by means of several techniques ranging 
from magnctorcsistivity 3 -^ to ac magnetometry itself 3 ^ 
and nuclear magnetic resonance^ This feature has al- 
ready been reported in the SmFeAsOo.8Fo.2 sample 1 ^ 
and, as it will be discussed later on, it can be justified 
in terms of pinning effects on a single FL propagating 
among bundles of entangled FLs. In this framework, the 
crossover between the two different trends of (Uo(H)) pw d 
vs. H shown in Fig. [5]can be interpreted as the transition 
from a basically single-flux line response at low H values 
to a collective response of vortices for H > 5000 Oc. The 
saturated low- if values for the depinning energy barriers 
are typically (Uo(H)) pu ,d ~ 10 4 K. Such high values are 
in agreement with what reported from magnetoresistiv- 
ity measurements in SmFeAsOo.85 (see Ref. lilil ) and in 
Bai_a;K x Fc2As2 single-crystals, even if in the latter case 
a much weaker ii-dependence was observed^ These re- 
sults are clearly indicative of strong intrinsic pinning of 
the vortex lines. Features of strong-pinning mechanisms 
have also been deduced in single-crystals of PrFeAsOi-^ 
and NdFeAsOi-^-F^, confirming that the observed be- 
haviour is intrinsic in REFeAsOi^^F^ materials^ 



It must be noticed that, since the variation of 1/T p vs. 
v m is very modest (see the vertical scale in the inset of 
Fig. [5|), (Uo(H)) pw d is not only determined at a fixed 
H but almost in isothermal conditions as well. By re- 
ferring to the inset of Fig. [5J in fact, one can recognize 
that T p (v m ) is varying on a range of some tenths of K 
degree, as already observed in the case of the Sm-based 
sample. 13 This is a great advantage of the ac susceptibil- 
ity technique if compared, for instance, to magnetoresis- 
tivity measurements where (Uo(H)) pw d is estimated from 
an activated-like fit to data over a range of several tens 
of K degrees i 30 ' 40 One can then notice that the average 
temperature T* characterizing the variation of T p over 
the considered v m range is intrinsically positioned over 
the irreversibility line, so that the energy barrier should 
more correctly be referred to as (Uo(H* ,T*)) pw d, where 
(H*,T*) are the points on the H - T phase diagram 
belonging to the irreversibility line. 
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IV. ANALYSIS OF THE RESULTS 

The observed behaviour can be explained by means 
of the phenomcnological GL theory 3 ^ and by referring 
to the model of single-vortex pinning by atomic impuri- 
ties (like, for instance, ionic substitutions like 2_ /F~ 
or 2 ~ vacancies) i 20 i 21 ' 37 Due to the very small values 
of the coherence lengths ~ 10 A, in fact, local defects 
on the atomic scale can be considered as strongly effi- 
cient pinning centers for FLs4^ In this framework the 
energy Uo is related to the FLs features only and does 
not depend on the precise pinning mechanisms.— More- 
over, at strong fi-values the high density of FLs gives 
rise to entangled bundles of vortices around the central 
one physically bound to the atomic defecti 20 i 21 ' 37 Accord- 
ingly, the characteristic energy Uq can be directly linked 
to the geometrical properties of the flux line lattice and, 
in particular, to the typical volume V of the correlated 
vortex lines. The following phenomcnological expression 
can be envisaged (where Uq is expressed in K ) 20 i 21,37 



Uo(T) 



H 2 C (T) 



V_ 

k B 



(9) 



where the term between curly brackets quantifies the In- 
dependence of the superconducting condensation energy 
density. Two limiting cases can be considered. 

On the one hand, by gradually increasing if the den- 
sity of vortices steadily increases and, accordingly, cor- 
relations among vortices increase too. A crossover to a 
regime where the depinning process leads to a collective 
response of an increasing number of vortices is then ex- 
pected. In the simplified scenario of a square Abrikosov 
lattice of vortices, the quantity d = yj $o/if estimates 
the mutual distance among nearest neighbouring FLs. 
One can then assume that the correlations roughly ex- 
tend over a cylindric volume whose radius is given by d. 
Namely one has V ~ 7r£(T)<i>o/if , the characteristic size 
along the third dimension being determined by the coher- 
ence length. 37 After the substitution in Eq. ©, consid- 
erations on the powder-averaging making it possible to 
obtain a more convenient form in order to describe the 
experimental data. In particular, the following relation 
holds under a two- fluid approximation (see Appcndix IA II 
for details) 



(Uq(H* ,T*)) pw d 



5/2 



g(t*)yJ(H c2 (0)) pwd 



1 

/i(7)A 2 b (0) 



if* 



(10) 



where g(t) is a function of t and /i (7) is a function of the 
anisotropy parameter 7 = S, a b/S,c = K/Kb- 

On the other hand, at low enough magnetic fields the 
typical volume V will be no longer sensitive to the pres- 
ence of several flux lines but it will only be a function 
of the typical lengths of the single pinned vortex. As a 
consequence, one can consider a different cylindric vol- 
ume with radius <5£(T). Here the heuristic parameter 5 is 



introduced in order to grant a continuous crossover be- 
tween the high- and the \ow-H regimes^ 7 - Thus, again 
choosing the coherence length as the characteristic size 
also along the third dimension, the following relation fol- 
lows 



(U (H*,T*)) pwd = 



h{t* 



192v^ 5 / 2 fc B y/(H c2 (0)) pwd 
1 

/ 2 (7)A 2 b (0)' 



(11) 



Here h(t) is a function of t and ^(7) is a function of the 
anisotropy parameter (details relative to the derivation 
of Eq. (TTTj) can be found in Appendix IA2[) . 

The two resulting expressions show that in both the H- 
regimes A a f,(0) can be simply expressed as a function of 
the experimentally-accessible quantity (Uq(H*, T*)) pwd . 
One should consider that the continuity of A a h(0) must 
clearly be assured at the crossover fieldi 3 - 7 - The free pa- 
rameters for the three different samples (namely, the 
anisotropy parameter "/re and the parameter 5 re), any- 
way, make the fulfilment of this requirement quite arbi- 
trary. Some other criterion from the analysis of experi- 
mental data should be formulated in order to derive some 
relations among the parameters, reducing as much as pos- 
sible the degree of arbitrariness. This can be done by the 
examination of the FLs phase diagram and by taking into 
consideration the role of thermal fluctuations. As a final 
result of the procedure, as reported in detail in Appendix 
[Bl one finds that 



3 / 2 x2 1 a 3 / 2 x2 o 3 / 2 x2 

7Sm 5 Sm ~ l- 6 7 Ce S Ce - ?>1 La 5 La- 



(12) 



The linking up the different sets of data can now be per- 
formed under the fulfilment of the constraints reported 
in Eq. (fT2|) . As a starting point of the procedure, a 
reasonable value for one of the "/re parameters should 
be taken as fixed. By referring to typical data reported 
in literature f* 3 - in particular, it will be assumed that 
Ism = 5. Next, data relative to the two different regimes 
(high- and low- if) in SmFeAsOo.sFo.2 are linked up by 
setting a proper value of Ss m - The procedure is repeated 
also for Ce- and La-based samples, where the starting 
guesses for the parameters "/re and Sre must be modi- 
fied till the fulfilment of Eq. (15). 



V. DISCUSSION 

As reported in the inset of Fig. [5J the linking pro- 
cedure of data in the two opposite limits described by 
Eqs. p7J|) and (TTJ) under the constraints reported in 
Eq. (TJJ allows one to deduce a sizeable if-dependence 
of X a b(0) itself. The saturation values in the limit 
if — > Oe have been reported in the last column of 
Tab. |TT] compared to the results of /x + SR measurements 
on the same samples (raw data not shown), together 
with the results for jre and Sre- The presented results 
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H(Oe) 

FIG. 6: (Color online) Main panel: graphical representation 
of the collapse of experimental data represented in Fig. [5] af- 
ter a proper normalization according to the value of 1/X ab (0) 
estimated at 250 Oe. The continuous line is a best-fit to data 
according to Eq. (|13[) . Inset: //"-dependence of the A a f,(0) val- 
ues for the three investigated samples. Open (full) symbols 
are relative to data from the low-H (high-//) regime. 



of /i + SR measurements are in good agreement with re- 
sults on different samples characterized by slightly differ- 
ent stoichiometrics 1 44 ' 45 At the same time, the estimates 
of X a b(0) for CeFeAsO0.92F0.08 and SmFeAsOo.sFo.2 by 
means of \ac are consistent with the values from /i + SR 
measurements within a systematic factor ~ 2. This dis- 
crepancy could possibly be accounted for by assuming 
that the correlations among vortices in the high-// regime 
are actually extended over a bigger volume. In particu- 
lar, it is immediate to realize that the choice of 2d for the 
cylindric volume, also implying a rescaling 5 — > 25, leads 
to a complete agreement of \ac data with ^ + SR ones. 

The observed behaviour of 8 as a function of the RE 
ion can be qualitatively understood in terms of a com- 
pensation of the opposite trend in (^) pw d- The trend in 
the modification of the value of 7, on the other hand, can 



be correlated with the increase in the size of the liquid 
region in the phase diagram reported in Fig. |4l One 
can deduce that the increase in 7 possibly leads to the 
enhancement of 2-dimensional fluctuations, much more 
effective than the higher-dimensional ones in extending 
the liquid region of the phase diagram. It should be re- 
marked that a similar trend in La- and Sm-based samples 
was reported in literature, even if the measured absolute 
values were considerably higher 4 s - 

More interestingly, a scaling of A~ h 2 (0) is displayed in 
the main panel of Fig. [6] showing a clearly common in- 
dependence after a proper normalization with respect to 
the value A~ fe 2 (0) taken at H = 250 Oc. This feature pos- 
sibly demonstrates the existence of a common material- 
independent underlying behaviour. It should also be re- 
marked that the quantity A~ fc 2 (0) is directly associated 
with the superfluid density n s of the superconducting 
state. 35 ' 36 This, in turn, implies that n s is partially sup- 
pressed by H values much lower than H C 2 (0) , a character- 
istic fingerprint of multi-gap superconductivity . 10 ' 47 By 
assuming a //-dependence typical of s-like bands^ one 
can derive the following phenomenological fitting func- 
tion 

n s {H) = n sl + n s2 {H) = n sl H „ . (13) 

The fitting results show that the main contribution to 
the superfluid density (1 — n s i) ~ 0.9 comes from the 
weakest band that is suppressed by a typical magnetic 
field Ho ~ 1.25 x 10 4 Oe. This quite unphysical result 
should possibly be associated with the high degree of ap- 
proximation associated with the two-fluid model and the 
relative T-dependence of the characteristic lengths (see 
Eqs. (|A4I) later in Appendix I A 1| . Our results are any- 
way in good qualitative agreement with previous reports 
obtained by means of torque magnctometry and local 
muon spin spectroscopy on an optimally-doped Sm-based 
superconductor.— In that case from ^t + SR it was possible 
to deduce a weaker //-dependence, eventually saturating 
at [IqH ~ 1.5 T and involving the suppression of just 
25% of the overall superfluid density. 



TABLE II: Parameters 5, 7 and X a b(0) (estrapolated for 
H — > Oe) for the three samples resulting from the con- 
strained linking procedure of data relative to the two distinct 
//-regimes (see text for details). Values of A a ;,(0) as obtained 
from /i + SR measurements have been reported in the last col- 
umn (raw data not shown in the text). 



Sample 


7 


S 


A a6 (0) (nm) 


A a6 (0) (nm) 








from Xac 


from p + SR 


LaFeAsOo.gFo.i 


3.2 


9.5 


235 ± 15 


/ 


CeFeAsO0.92F0.08 


4.25 


10.5 


128 ± 6 


260 ± 15 


SmFeAsOo.sFo.2 


5 


11.75 


98 ± 8 


200 ± 15 



VI. CONCLUSIONS 

The phase diagram of the flux lines in three powder 
samples of REFcAsOi-zF^ superconductors (RE = La, 
Ce, Sm) under conditions of nearly-optimal doping was 
investigated by means of ac susceptibility measurements. 
The irreversibility line has been estimated for the three 
samples, showing that in the accessible range of mag- 
netic field the extension of the reversible liquid region is 
lowered with increasing T c . This aspect could be of ex- 
treme interest in view of possible technological applica- 
tions of these materials. The //-dependence of the char- 
acteristic depinning energy barriers was investigated in a 
thermally-activated framework. Results were interpreted 



8 



by distinguishing two regimes of magnetic field character- 
ized by single and collective depinning processes, allowing 
us to make experimental data collapse on the same tem- 
perature trend indicative of a common underlying mech- 
anism independent on the precise material. Reliable esti- 
mates of the zero-temperature penetration depth A a b(0) 
and of the RE-dcpcndcncc of the anisotropy parameter 
7 are finally given. Some interpretations of the observed 
behaviour in terms of a two-band s^ 1 model are proposed 
even if a refinement of the employed two-fluid model is 
in order. 



Acknowledgements 

R. Khasanov, G. Lamura and A. Rigamonti are grate- 
fully acknowledged for stimulating discussions. H. Stum- 
mer, C. Malbrich, S. Muller-Litvanyi, R. Muller, K. 
Leger, J. Werner and S. Pichl are acknowledged for assis- 
tance in the sample preparation. The work at the IFW 
Dresden was supported by the Deutsche Forschungsge- 
meinschaft through Grant No. BE1749/12 the Priority 
Program SPP1458 (Grant No BE1749/13). S. W. ac- 
knowledges support by DFG under the Emmy-Noether 
program (Grant No. WU595/3-1). 



Appendix A: Derivation of the relations between 

{Uo(H* , T*)) pw d and A a &(0) in the two different 
H-regimes 

1. Strongly correlated vortices. Depinning of 
bundles of flux lines 

Let's consider Eq. © in the limit of high-i/, here 
rewritten for convenience 



U (H,T) = 



Hl{T)j{T)^ 
8k B 



H 



(Al) 



as the starting point. By means of the GL relation for the 
flux quantum $ = 2\/2ttH c (T)\(T)£(T), one can write 
down an explicit expression for Uq in different conditions 
of orientation of the magnetic field as follows 



U ,H\\ab(H,T) = 

U . mc (H,T) = 



1 



^ 2 k B H y/U(T)S c (T)\ ab (T)\ c (T) 

(A2) 



&A^k B H U(T)\ 2 ab (T)- 



The expressions £ a b/£ c = A c /A a 5 = 7 hold for the differ- 
ent typical lengths in anisotropic superconductors, where 
7 has already been defined in the text as the anisotropy 
parameter. By referring to Eq. ([1} it is possible to per- 
form a powder-like average of Uq as 

(U (H,T)) pwd = lu , m[ab (T) + ±U 0Mc (T) = (A3) 



d>3 



1 



V7 + 2 

3 7 3/2 J U^k B H UT)K b {TY 



In a simple two- fluid model, the T-dcpcndcnce of £ C (T) 
and \ a b(T) can be taken as 



fcCn = &(0)^O-, A a6 (T) = A o6 (0)- 7 =L= (A4) 

where t has already been defined as the reduced tem- 
perature t = T/T c (0). It should be considered that, as 
already recalled, the estimate of (Uo(H,T)) pw d is per- 
formed by definition at the values (H*,T*) delimiting 
the irreversibility line. After the definition of the func- 
tion g{t) = (l — i 2 ) y/1 — t 4 it is then possible to write 



(U (H*,T*)) 

pwd 

g(t*) 



<r>ii 



V7 + 2 

3 7 3/2 J 647r 2 fc B -H"* 
1 



(A5) 



The quantity £ c (0) can be independently derived by mea- 
suring the magnetic field dependence of the supercon- 
ducting transition temperature T c by means of the rela- 
tions reported in Eqs. (|S|) and ([5]) after proper consid- 
erations about the powder-average procedures. The two 
limiting configurations H \\ ab and H \\ c lead to the 
formulas 



H c 2,H\\ab{T) 
H c 2,H\\c(T) 



27T&CT)' 



(A6) 



Following Eq. ([T]) it is possible to deduce that the 
experimentally-accessible quantity (£{0))pwd) already de- 
fined in Eq. ((6]), is linked to £ c (0) by the relation 



Coming back to Eq. (|A5[) . one can substitute Eqs. 
and ([5]) to obtain 



(A7) 



pwd 



<i> 



5/2 



g(t*)^/(H c2 (0)) pwd 



96\/27r 3 / 2 fc B 
1 1 



" /i(7)A 2 ,(0)^ (A8) 
having defined the function of the anisotropy parameter 



Mi) = 



1 + 7 



Eq. 



V3 + 2 
is equivalent to Eq. (|TU|) . 



(A9) 



2. Weakly correlated vortices. Depinning of single 
flux lines 



Let's now consider Eq. (JSJ in the limit of low- H, here 
rewritten for convenience 



Uo(T) 



(AlO) 
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Again by means of the GL relation for the flux quantum 
$o = 2\/27r H C (T)X(T)£(T) it is possible to explicit the 
expressions for Uo in the different cases of the orientation 
of the system with respect to the magnetic field as 



Uo,H\\ab(T) 
U ,H\\c(T) 



64n 2 k B X ab (T)X c (T) 

<^ 2 U{T) 
ttePk B x 2 ab (Ty 



(All) 



Similarly to what performed in the previous Appendix 
concerning strongly correlated vortices, the 7 factor can 
be introduced and X a b(T), £ C {T) can be left as inde- 
pendent quantities. By considering that the estimate of 
(Uo(H 7 T)) pw d is performed along the irreversibility line 
(H*,T*), by employing Eq. (|A4[) and after a powder-like 
average of the energy barrier one can write 



(Uq(H*, T*)) pw d 

W) 



7 2 + 2/7 
6(0) 

KM 



64ir 2 k B 



(A12) 



where h(i) = (l + t 2 ) yl — t A . By again expressing £ c (0) 
in terms of (H C 2(0)) pw d through Eq. (|A7|) and by means 
of Eq. © , it is finally possible to deduce the following 
expression 



(Uq(H*, T*)) pw d\J (H 02(0)) P wd 



192^7T 5 / 2 fc B 

1 



(A13) 



(A14) 



In the previous expression, f 2 {^) is defined as 

/2(7) = 

32K1) ( 7 2 + 2^7) 

Eq. (|A13|) is equivalent to Eq. (fTT|) . 



Appendix B: Linking procedure of data in the two 
different //-regimes 

In this Appendix the problem of the continuity at the 
crossover field H cr of the A a h(0) vs. H data, obtained 
by means of Eqs. (fTU)) and (fTTj) , will be considered. In 
particular, as already stated in the text, some constraints 
on the variability of the six parameters jre and <5 re (RE 
= La, Cc, Sm) should be fixed in order to reduce as much 
as possible the degree of arbitrariness of the procedure of 
data-linking. 

At this aim, it is convenient to introduce the Ginzburg- 
Levanyuk number Gi 3 o{0) as^ 



2 fc B T c (0)A 2 (0) 



(Bl) 



quantifying the extension of the region of the H — T 
phase diagram where thermal fluctuations are sizeable 
and significantly affect the physics of the system. The 
//-dependence of Gi 3 o(0) is given byi£ 



Gi 3D (0,H) = 



2H^Gi 3D {Q) 
H c2 (0) 



2/3 



(B2) 



In fact, one can assume that the position of the ir- 
reversibility line is mainly governed by the amount of 
thermal fluctuations in the system. As a consequence, 
Gi3u(0, H) is expected to be directly involved in the an- 
alytic expression relative to the irreversibility line itself. 
One hint at the correctness of this picture is possibly 
given by the similarity between the characteristic expo- 
nents observed in Eqs. ((7|) and (|B2j) . Similar consider- 
ations, moreover, have already been proposed in litera- 
ture concerning the thermodynamical melting line (see, 
in particular, Sections IV and V of Ref. l5(iland references 
therein. In that case, anyway, the considered exponent 
is j3 = 2). The following expression for the irreversibility 
line can then be considered 



Tc(0) 



«3D(0,g*) 

K 2/3 



(B3) 



where K is an arbitrary proportionality factor. Together 
with Eq. (|B2[) . this straightforwardly leads to 



2H*^Gi 3D {Q) 
H c2 (0) 



= K 



1 - 



T c (0) 



3/2 



(B4) 



Since all the sample-dependent quantities are already 
kept into consideration by Gi 3 d(0)-, it is reasonable to 
assume K as a sample-independent parameter. K can 
typically be interpreted in terms of microscopic proper- 
ties of the vortices £2- In the present phenomenological 
model, anyway, such microscopic interpretations are left 
aside. 

In order to give a suitable description of the exper- 
imental data, Eq. (|B4[) should be powder-averaged by 
considering the criterion already presented in Eq. (J]). 
By considering the definition of Gi 3 d(0) reported in Eq. 
(|Bip , one can write 



Gi 3D (0) 



1 8^ 2 fc s T c (0)A ab (0)A c (0) 27r£ a6 (0)£ c (0) 



1 8^ 2 fc B T c (0)A 2 b (0) 2< 2 fc (0) 



Gi 3D (0) 
H c2 (0) 



V2 *§U(0) 



(B5) 



for the two field orientations H \\ ab and H \\ c, re- 
spectively. By again performing the powder-average pro- 
cedure and after introducing the anisotropy parameter 
through the relations S, a b/^c = X c /X ab = 7, by referring 
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FIG. 7: (Color online) Graphical representation of the quan- 
tity between square brackets in Eq. (|B9I) . The dashed lines 
are guides for the eye and they can be used to extrapolate the 
intercept values back to H = Oe. 



to Eqs. (|6|) and (|A7|) one obtains the following expression 

1 8i: 2 k B T c (0)X 2 ab (0) 



Gi 3D (0) \ 
c2 (0) / 



H, 



pwd 



Hi) V2*i&(0) 
i 



(H c 2(0))pwd 

where the function 73(7) is defined as 

9 7 2 



/s(7) 



(2 7 3/2 + 7 ) (2 7 + 1) ' 



(B6) 



(B7) 



One should now consider that the definition of GisoiO) 
reported in Eq. (|B1[) docs not account for any in- 
dependence that is fully accounted for by Eq. (|B2j) . Eq. 
(|B6j) is then referred to a H = Oc condition and, as a 
consequence, the quantity A 2 fc (0)/£ c (0) can be obtained 
by means of Eq. (|A12[) (holding in the low-H regime). 
By inserting the resulting expression into Eq. (|B4j) one 
obtains 



H* 



\H c 2(0))pwd 



K RE X 



0.95 



Tc(0) 



3/2 



where 



K RE = 



_ / (U (H*,T*)) pwd 



T c (0)h(F) 



27^2 



RE 



3/2 c9 

Ire 6 * 



(B8) 



(B9) 



RE 



As already discussed in the text, the factor 0.95 in Eq. 
(|B8p phenomenologically accounts for the definition of 
the irreversibility line from the Xa' c -critcrion. 

The quantity between curly brackets in Eq. (|B9|) can 
be experimentally estimated from a linear extrapolation 
of the intercept values back to H = Oe (see Fig. [7]). 
One can then compare the sample-dependent quantity 
Kre with the corresponding experimental quantities de- 
rived from the fitting procedure shown in Fig. [4] To- 
gether with the already cited assumption that K is a 
sample-independent quantity, this implies that some con- 
straints can be put on the variability of the parameters 
jre and 5 re- In particular, one finds that 



3/2 ,-2 
as reported in Eq. 



I.67, 



3/2,2 

Ce °Ce 
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